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ABSTRACT 



Context. Stellar jets from young stars can be highly asymmetric and have multiple velocity components. 

Aims. To clarify the origin of jet asymmetries and constrain the launch mechanism we study, as a test case, the physical and kinemat- 
ical structure of the prototypical asymmetric flow emitted by DG Tau B. 

Methods. The analysis of deep, high spectral resolution observations taken at the KECK telescope allows us to infer the properties and 
the spatial distribution of the velocity components in the two jet lobes. From selected line ratios we derive the gas physical conditions 
(the electron and total density, n e and n H , the ionisation fraction, x e , and the temperature, T f ), as a function of both distance from the 
source and gas velocity. The presence of dust grains in the jet is investigated by estimating the gas-phase abundance of calcium with 
respect to its solar value. 

Results. The detected lines show broad velocity profiles at the base of the jet (up to ~100kms _1 ) where up to three velocity compo- 
nents are detected. At 5" from the source, however, only the denser and more excited high velocity components survive and the lines 
are narrower (~ 10-30 kms -1 ). The jet is strongly asymmetric both in velocity and in its physical structure. The red lobe, which is 
slower (~140kms~') and more collimated (opening angle: ff~3-4°), presents low ionisation fractions (x e ~0.1-0.4) and temperatures 
(T e <5 10 3 K), while the total density is up to ~2.5 10 4 cm 5 . The blue lobe, faster ( — 320 kms~') and less collimated (a~l4°), is 
also less dense (n H ~l 10 4 ever 3 ) but highly excited (T e up to ~5 10 4 K and x e up to 0.9). The estimated mass loss rate turns out 
to be similar in the two lobes (~6-8 10~ 9 M Q yr~' ), while the flux of linear momentum is 3 times higher in the blue one (~2.5 10~ 7 
M G yr" 1 km s~ x ). Calcium is strongly depleted with respect to its solar abundance, indicating that the jet contains dust grains. The 
depletion is lower for higher velocities, consistent with dust destruction by shocks. 

Conclusions. The similar mass loss rate in the two lobes suggests that the ejection power is comparable on the two sides of the system, 
as expected from a magneto-centrifugal ejection mechanism, and that the observed asymmetries are due to different mass load and 
propagation properties in an inhomogeneous environment. The presence of dust grains implies that the jet is generated from a region 
of the disk extending beyond the dust sublimation radius. 

Key words. ISM: jets and outflows - Herbig-Haro objects - dust, extinction - Stars: formation 
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1. Introduction 

Stellar jets from young stellar objects (YSO) play a key role 
in the star formation process and a wealth of information can 
be derived analysing their characteristic emission line spectra. 
The observed lines are collisionally excited in the shock waves 
generated by the interaction of the jet material with the inter- 
stellar medium or previously ejected matter and contain im- 
portant information on the gas physics/kinematics. Different 
methods have been proposed to derive the physical condi- 
tions of the gas propagating in the jets from emission line ra- 
tios, such as the compar ison with the predic t ions of a grid 
of sh ock models (e.g., |Raga & Bohm 1986; Hartiga n~et alJ 



1994. or, alternatively, sp ectral diagnostics techniques (e.g., 
Bacciotti & Eisloffel 1999, hereafter referred to as the BE tech- 
nique). 



The application of these methods to analyse emission lines 
allows one to derive the jet structure on different scales, 
depending on the angular resolution of the observations : 
from parsec scales to hun d reds of AU jHartigan et al. 1994; 
Bacciotti & Eisloffel fl999t iPodio et alJ 120061) . and down to 
~15 AU from the emitting s ource with space or A d aptive 
Optics-assisted obse r vations jLavallev-Fouauet et ail 2000; 



2007; Melnikov et al 



Bacciotti et al. 2000, Woit as et all 120021; lHartigan & Morsel 



2009). In particular, the analysis of emis- 



sion lines can shed light on the gas conditions at the base of 
the flow and on the mechanism that generates it. Proposed mag- 
neto hydro-dynamical (MHD) models suggest that the wind is 
accelerated and collimated by magneto-centrifugal forces. It is 
not clear, however, where the jet originates: from the star it- 
self (iSautv & Tsin ganoslll994l) . from the radius at which the 
stellar magnetosphere truncates the disk (X wind, IShu et alJ 
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2000), or from an ext ended region of the disk (Disk wind, 



Konigl & Pud"ritzll2000l) 



In addition to the uncertainty on the launch mechanism, 
other observed properties still lack a proper modeling. For ex- 
ample, velocity resolved observations of the jets from T Tauri 
stars and younger Class 0/1 protostars show that these can 
be hi ghly asymmetric and can have different velocity compo- 
nents dHamandll994lffirth et al.ll 1 994 1 1 997h iDavis et alJl200lt 
iGarcia Lopez et al.ll2008l 1201 Oh . Despite some attempts to ex- 
plain these ob servational features by means of the exis ting MHD 
models (e.g., iPesenti et al.1 12003: Ferr eira et al.l l2006) their ori- 
gin is still unclear. A viable possibility to clarify these as- 
pects is to investigate the variation of the jet physical proper- 
ties in asymmetric jets and in different velocity channels, as it 
has been attempted in recent wor ks conducted at moderate and 
high spatial/spectral resolution ([L avallev-Fouquet et al. 2000; 
Woit as et alJl2002t ICoffev et al.ll2008t IGarcia Lopez et al]l2008L 
2010; lMelnikov et al.ll2009HPodio et alj|2009l) . 



Another interesting issue, which remains not fully clarified, 
is the dust content of jets. Refractory species such as Ca, Ni, 
Cr, Fe, are often locked onto dust grains, thus a strong de- 
pletion of their gas ph ase abundance is expecte d in the inter- 
stellar medium (ISM) (ISavage & Sembachlll996l) . A few stud- 
ies have i nvestigated the gas phase abundance of Fe and Ni in 
HH j ets (iBeck-Winchatz et all 1 1994 Il 996t iMouri & Taniguchl 
2000; Bo hm & Matlll200lt iNisini et al.ll2002l) and only recently 
the analysis has been ext ended to other refractory species such 
as Si Ca, C, Cr, and Ti dNisini et al.1 120051 l2007t ICabrit et all 
120071 lPodio eTaT] l2006. 20091: IGarcia Lopez et al.ll2008l) . These 
studies showed that the refractory species may be depleted up 
to 90% in HH jets and molecular outflows, thus suggesting the 
presence of dust grains in the ejected material. This in turn can 
put constraints on the region of the disk from where jets origi- 
nate. In fact the stellar radiation destroys dust grains in the disk 
up to the so called dust evaporation radius, R evp . This is located 
between 0.03 and 1.5 AU from the star in T Tauri and Herbig 
Ae/Be star s depending on the stellar luminosity and the dust 
properties JlseTl a & Natta fl2005tlAkeson et al.l2005llEisner etall 
I2007T) . but no direct observations and/or modeling of the inner 
disk structure are available for younger Class 0/1 sources, still 
surrounded by circumstellar matter. In general, we expect that 
a 'dusty' jet can only arise from a region in the disk extended 
beyond R evp . The passage through a shock front, however, can 
destroy part of the dust grains wh ich may be transported in the 
jet, co mplicating the picture (e.g., May et al. 2000; Guille t et al.l 
2009). Therefore, accurate studies of the gas-phase abundance 
of refractory species are necessary to determine the dust content 
in the jet and to constrain the shock efficiency in destroying the 
dust grains. 

To improve our understanding of the physics of HH jets from 
YSO, particularly their launch mechanism, we investigate the 
physical and kinematical structure of the HH 159 jet emitted by 
the young source DG Tau B, focusing on the analysis of the ve- 
locity components detected at its base, the asymmetry between 
the two jet lobes, and, finally, the estimate of the dust content in 
the jet. 

HH 159 is a bipolar jet, first detected by Mun dt & Fried! 
(119831) . The red lobe consists of a chain of bright knots extend- 
ing to ~55" from the source, while the blue lobe, fainter and 
less collimated, is detected only up to ~10 " from the source 
dMundt et al.lll99lHEisloffel & Mundtlll998l) . A large molecular 
outflow spatially coincident with the redshifted optical jet ha s 
been detected in the CO lines by iMitchell etaTI d!994 1 19971) . 
The driving star, DG Tau B, has been classified as a Class I 



source based on its spectr al energy distribution dWatson et al.l 
2004: iLuhman et aT1l2010l) and at optical wavelengths it is ob- 
scured by circumstellar optically thick material detected in ab- 
sorption through broad band imaging wi th the Hubble Space 
Telescope (HST) dStapelfeldt et all 1 1997b . The observed dark 
lane, elongated and perpendicular to the jet, indicates the pres- 
ence of a c ircumstellar disk, confirm ed by 13 CO millimeter ob- 
servations dPadgett et al.ll 19991 120001) . and flattened residual en- 
velope. The source position has been determined from 3.5 cm 
VLA radio continuum observations by Rodrigu ez et al.l d 19951) . 
who locate it within the observed dark lane. Further HST studies 
in the near-infrared indicated the presence of a bipolar reflection 
nebula, whose axis of symmetry is coincident with the axis of 
the optical jet (Pa dgett et al.| [T999). The eastern lobe of the neb- 
ula, surrounding the blue jet lobe, is V-shaped, suggesting that 
the nebula traces the walls of the blueshifted outflow cavity. The 
western lobe is fainter and more collimated, and encompasses 

the redshifted optical jet and the CO ou tflow. 

dMundt et al 



Previous 



1987 



1991 



s spectroscopic 
Eisl offel & Mundtl Il998l) and imaging dMundt et al.l 
Stapelfel dtet alJ Il997t iPadgett et al. 1999) studies highlighted 
a strong asymmetry between the two lobes in morphology, 
velocity and degree of excitation. These studies, however, 
focused only on the Ha and [S n] emission lines as well as on 
the continuum emission in the optical and near-infrared bands. 
On the contrary, in this paper we present for the first time very 
deep and high-spectral resolution spectroscopic observations of 
a large number of forbidden and permitted lines in the optical 
range from 5300 A to 8500 A. The high spectral resolution 
coupled with long integration times highlight a very complex 
velocity structure with multiple velocity components evolving 
along the jet. The large number of detected lines has allowed 
a detailed study of the gas physical conditions, through the 
application of spectral diagnostic techniques to selected line 
ratios. This has provided a rich information on the excitation 
and dynamics of the flow, as well as on the dust reprocessing. 
The comparison of our results with those obtained for other 
jets provides means of investigating the nature of the detected 
velocity components and the origin of asymmetric jets. This, 
in turn, can help understanding the generation of outflows and 
their role in young systems. 



2. Observations and data reduction 

High resolution echelle spectra of the jet from DG Tau B were 
obtained with the HIRES spectrograph at the Keck I telescope 
dVogt et al.lll994l) . located at the Keck Observatory on Mauna 
Kea, Hawaii, on February 23/24, 1998. The weather conditions 
were stable and photometric during the nights, providing a see- 
ing of typically 077-078. Four deep spectra of 2700 s integration 
time were obtained using a long slit decker of 28"xl715, and a 
2048x2048 CCD Tektronix detector with 2xbinning in the spa- 
tial direction. This setting provided very high spectral resolution 
(R~34 000, i.e. -8.8 km s" 1 ), spectral sampling of ~0.036-0.052 
A pixel, (corresponding to ~2.15 kms -1 ), and a spatial scale 
of -07378 pixeF 1 . 

The d ecker was aligned alon g the jet position angle 
(PA~295°, lEisloffel& Mundtl 119981). Since the deck er is 28" 
long, i.e. shorter than the jet (jEisloffel & Mundt 1998), to cover 
the full jet length the observations were carried out by centering 
the decker on the jet source, and then shifting it along the jet 
PA towards the south-east and north-west relative to the source 
position, along the blue and the red lobe. In this way a length 
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Fig. 1. Position- Velocity diagrams of the jet from DG Tau B in the strong optical emission lines [S n]/l6731, [Oi]/I6300, [Nn]/16583, 
and Havl6563 (angular resolution: FWHM, f , e ,„ ( ,~0.8", spectral resolution: FWHM, nj ,~8.8 kms -1 ). Contours are equally spaced in 
logarithmic scale. The first contour is at 3cr level for the [S n], [Oi], and [Nn] lines and at 5<x for the Ha line. The peak level is 
at 1 103cr for the [S n], 332cr for the [Oi], 1 15<x for the [Nn], and 131<x for the Ha line. The vertical and horizontal dotted lines 
indicate the zero radial velocity with respect to the star and the spatial position of the source. The dashed horizontal lines indicate 
the position of the detected knots. 



of ~55" was covered along the spatial direction. The position 
angle tracking mode allowed us to keep the slit orientation at 
the chosen PA during the four exposures. By using the decker, 
23 spectral orders overlapping along the spatial direction were 
obtained, covering a wavelength range from 5300 A to 8500 A. 

The raw frames were reduced following standard recipes 
with the noao.imred.ccdred package in IRABJ, including over- 
scan correction, bias and dark subtraction, flat-fielding, back- 
ground subtraction, and wavelength calibration using Thorium- 
Argon arcs. The Laplacian Cosmi c Ray Identification method 
(iacos' method. Ivan Dokkumll200l1) was used to clean the im- 
ages from cosmics. 



1 IRAF is distributed by National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc. under contract with the National Science Foundation 



Because of the characteristic jet emission line spectrum we 
are interested in short sections of the echelle spectra. Thus, the 
latter were transformed into ordinary long-slit spectra centered 
around the detected lines. A flat-field frame obtained with a 28" 
long slit was used to calculate the tilt of the orders in the echelle 
images with respect to a normal long slit spectrum. Hence, the 
orders containing the emission lines were corrected for the es- 
timated tilt. The same procedure was applied to the Thorium- 
Argon arcs frames. 

Wavelength calibration for the selected spectral orders was 
performed by means of Thorium- Argon arcs, providing an accu- 
racy of about 0.003 A. The observed velocities were corrected 
for the velocity of DG Tau B with respect to the observer, hence 
the values of V m d showed in all the plots refer to the gas ra- 
dial velocity with respect to the driving source. Since the central 
source is obscured it is not possible to estimate the source ra- 
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Fig.2. Position- Velocity diagrams of the jet from DG Tau B in the [Ni],U5198, 5200, [Fen]^5159, [Nin]^7378, [Can]^7291 
emission lines. Note that the velocity scale in the first panel is relative to the [Ni]/i5200 line, while the scale for the [Ni]/L15198 
plotted in the same panel is shifted by -136 kms -1 in velocity. Contours are equally spaced in logarithmic scale. The first contour 
is at 2<x level. The peak level is at 37<x for the [N 1], 24<x for the [Fe 11], 37cr for the [Ni n], and 80<r for the [Ca n] line. The vertical 
and horizontal dotted lines indicate the zero radial velocity with respect to the star and the spatial position of the source. The dashed 
horizontal lines indicate the position of the detected knots. 



dial velocity fro m photospheric absorption lines, e.g. Li 6710 A. 
Thus, following Eisloffel & Mundtl |l998) we assume that DG 
Tau B has the same heliocentric ra dial velocity as the su rround- 
ing molecular gas (~16.3 kms -1 , Mitchel l et al] 1 19941) . which 
translates in a source velocity with respect to the observer of 
~46.4kms _1 . 

Finally, the four frames centered on different positions along 
the jet PA were aligned and then combined by estimating the po- 
sition of the brightest emission knots. Thus we obtain for each 
detected line a long slit spectrum covering ~55" along the spa- 
tial direction and ~3 10 kms 1 along the spectral one (see the 
position-velocity (PV) diagrams in Fig. Q] and [2] for selected 
lines). Note that in all the figures shown in this paper the dis- 
tance from the source is positive along the blue lobe and negative 
along the red one. 

At the end of the data reduction process, telluric features 
overlapping with the jet emission lines and background caused 
by the overlap of adjacent orders were fitted by a polynomial 
surface and then subtracted. The extremely bright [O i]/l6300 tel- 
luric line could not be properly subtracted, thus the [Oi]/16300 
line PV shows a ~20kms _1 cut in the blue lobe centered at — 
46.4 km s -1 (see Fig.Q]). This cut is also clearly visible in the PV 
diagrams of the derived line ratios and of the jet physical prop- 
erties (see Fig.[5]and Fig. [6). 

Note that the spectra are not flux calibrated, since no stan- 
dard star observations are available. Thus, the intensity of the 
lines in Fig.[T]and|2]is in arbitrary units with respect to the back- 
ground noise, cr. 

3. Results 

The deep KECK/HIRES spectra (see Fig.[T]and Fig. [2j allowed 
us to analyse the emission of the HH 159 jet up to ~40" in 



the bright red lobe ( knots Al, A3, A2, B2, C, and D, following 
the nomenclature bylHsloffel & Mundtll 1998b and -10" in the 
fainter blue lobe (knot Fl, F2, and G). Besides the [S \i\AA61\6, 
67 3 1 and Ha /16563 lines, wh ich have been previous ly observed 
bv lMundt et all d 19871 [19911) and lEisloffel & Mundtl Il998h . we 
detected strong and broad emission in [Oi]/t/l6300, 6364, and in 
[Nn]/l/i6548, 6583 (see Fig.[T]i. The jet is emitting also in a num- 
ber of fainter and narrower emission lines, such as [Fen]/I5159, 
[Ni]/U5198, 5200, [Caii]/I7291, and [Ni 11] ,17378, detected up 
to ~20" (knot C) (see Fig.|2]l. Finally, very faint emission in the 
Nai D lines at /L15890, 5896 is also detected at a ~6cr level in 
knots Al, A3. Since the ionisation energy of Nai is very low, 
i.e. ~5.1 eV, the Nai D lines are rarely obse rved in HH objects 
(see, e.g., HH 1 spectra by ISolf et alj[l988l) . The detection of 
such lines suggests that the excitation conditions are quite low 
in the red lobe. This idea is confirmed by our analysis of the gas 
physical conditions (see Sect. 13.21 ). 



3.1. Jet kinematics: asymmetries and velocity components 

The PV diagrams of selected emission lines in Fig.Q]show very 
broad (up to -200 km s FWZI) and complex (up to 3 velocity 
components) velocity profiles in the first few arcseconds from 
the source. Thanks to the high spectral resolution of the data 
(~8.8 kms -1 ) we can resolve and analyse each velocity compo- 
nent in the jet separately. To derive the kinematical properties 
we extracted the jet spectral profile in the two lobes by spatially 
integrating the [Sn]/I6731 and Ha emission lines over the see- 
ing full width at half maximum (FWHM Je , e ,„g~0"7-0'.'8, i.e. ~2 
pixels) at different distances along the jet length. Fig. [3] show 
the variation of the [Sn]/I6731 and Ha velocity profile in the 
red and the blue lobe up to ~6" from the source. The red lobe 
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Fig. 3. Normalized [S ii]/1673 1 (left panel) and Ha/16563 (right panel) velocity profiles along the first 6" of the blue and the red lobe 
(in arbitrary units), obtained from the line PV diagrams in Fig. Q] by spatially integrating the line spectral profiles over the seeing 
FWHM (~0.8", i.e. ~2 pixels). The label at each profile indicates the distance from the star. The vertical dotted line indicates the 
zero radial velocity with respect to the star. The dotted and dashed thick lines indicate the velocity peak of the high and the low (and 
moderate) velocity components, respectively. Note that beyond ~5", only the high velocity component survives in both lobes. 



clearly shows the presence of a low and a high velocity compo- 
nent (LVC and HVC), while in the blue lobe we detect three ve- 
locity components, at low, moderate, and high velocities (LVC, 
MVC, and HVC, respectively). Interestingly, in both lobes, the 
low and moderate velocity components fade away at ~4-5" from 
the source, while the high velocity components survive further 
out. This trend has been observed on a large sample of T Tauri 



stars but on much smaller scal es, i.e between a fe w te ns and 200 
AU fr om the source, by, e.g.. iHirth et alj ( 1 1997b and iPvo et al.l 
(120031) . and on similar or even larger scales at the base of HH 
jets from Class 0/1 sources dGarcia Lopez et alJ2008ll2010h . 

By deblending and Gaussian fitting the spectral profiles 
shown in Fig. [3] we obtain the gas radial velocity with respect 
to DG Tau B, V m &, and the spectral width of each velocity 
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Fig. 4. Jet kinematical properties as a function of the distance from the source. From top to bottom: (a) intensity profiles of the [S n] 
and the Ha lines integrated over their velocity profile and normalized to their peak intensity (red solid and green dotted-dashed line, 
respectively); (b, c) peak radial velocity, V ra d, and full width at half maximum, AV, of the different velocity components in the jet 
as inferred from Gaussian deblending of the [S ii]/16731 (HVC: red points, MVC, LVC: black points) and the He 6563 (HVCgreen 
points) line profiles shown in Fig. [3] 



component in the jet, AV, measured as the observed line full 
width at half maximum (FWHM i s ) corrected for the instrumen- 
tal width (FWHM,„ S ,~8.8 kms -1 ). This is done by assuming a 
Gaussian profile for both the instrumental and the emission-line 
profile (i.e. AV = (FWHNU., 2 - FWHM,„ sf 2 ) 1/2 ). The error on 
the obtained determination of V m d and AV varies between 1-4 
km s" 1 . The parameters derived by fitting the [S ii]/1673 1 and the 
Ha/16563 line profiles are plotted in Fig.|4]as a function of dis- 
tance from the source. 

The [Sn]/16731 line profiles in Fig. [3] show very well de- 
fined velocity components. In the red lobe the HVC radial 
velocity is decreasing in the first 14" from the source from 
~60kms~ I in the innermost knots Al, A3 down to ~47 kms 1 
in knots B2. Then V ra( j is increasing again moving towards the 
outer knots C (Vrad-SO-SSkms^ 1 ) and D (^-55-60^8^ ). 
This trend, clearly visible also in the line PV diagrams (see 
Fig. Q] and [2j, can be explained if the jet direction is chang- 
ing due to precession or interaction with the ambient medium. 
Observational evid ences of prec e ssing j ets have been alread y 
repo rted by, e.g., |G iannini et al. (2003J), iHodapp et alj (|2005), 
and ICaratti o Garatti et al.l (120081) . The LVC peak velocity in- 
stead varies between 25 and 42 kms _I . In the blue lobe much 
higher radial velocities are measured of -120-135 kms~' in 
the HVC, -20-30 km s" 1 in the MVC, and a few km s" 1 in the 
LVC. Given the jet inclination angle with respect to the line of 
sight of ~65° estimated t hrough proper motions measurements 
(Eisl6ffel&Mund3 ll998l) . the bulk of the gas is moving with 



a velocity Vj et ~140kms 1 in the red lobe, and more than two 
times higher (Vj e t~-320km s~') in the blue one. These va l ues ar e 
in agreement with those estimated by Eisl6f feT& Mundd (Tl998). 
The width, AV, of the high velocity component is decreasing 
with distance from the source in both lobes. In the blue lobe we 
measure AV~60 kms -1 at decreasing to ~30 kms _I at 
~6". In the red lobe AV decreases more steeply from ~40-60 
kms~' within 1" from the source down to -15 kms -1 at ~3", 
and values between 6-20 km s _1 for larger distances. 

The Ha /16563 line profiles are much broader than the [S n] 
profiles (FWZI up to ~200 kms _I ) and the emission is very dif- 
fuse in this line, making it difficult to separate the different ve- 
locity components. Thus we only report the values of V, a d and 
AV estimated for the HVC of the red lobe. Fig. [4] shows that 
V ra d nas the same trend as in the [S n] lines, with only slightly 
lower values. The line width is much larger, though (AV varies 
between ~35 and ~100 kms" 1 ). A bro ader Ha line profile was 
also noted bv lEisloffel & Mundtl (Il998l) . 

The large velocity dispersion measured within a few arcsec- 
onds from the source and, in particular, in the blue lobe, suggests 
that stronger shocks are present at the jet base, and in the blue 
lobe. As we will show in the next section this reflects in higher 
excitation conditions in the blue lobe, that is the jet is highly 
asymmetric not only in velocity but also in its physical structure. 
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Table 1. Kinematical, physical, and dynamical properties of the jet from DG Tau B in the high and the low velocity components 
(HVC and LVC, respectively) averaged over the first 5" of the red and the blue lobe. The values of V T& & and AV are derived from 
Gaussian fits to the [S n]/16731 line profiles and to the Ha 6563 (in parentheses). 



Lobe 


VC 


Vrad 


AV 


n e 


x f 




riH 


M jet 


-< jet 






kms~' 


kms~' 


10 3 cnT 3 




10 4 K 


10 4 cnT 3 


10-' M Q yr 1 


10- 6 Moyr'kms- 1 


Red 


HVC 


58 (54) 


27 (46) 


3.3 


0.22 


0.45 


1.49 


4.2 


0.6 




LVC 


35 


26 


1.9 


0.17 


0.40 


1.29 


2.2 


0.2 


Blue 


HVC 


-126 


47 


3.6 


0.88 


2.46 


0.41 


8.2 


2.5 




LVC 


-2 


18 


1.1 


0.41 


0.74 


0.29 


0.09 
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Fig.5. Position-velocity diagrams of the line ratios used to infer the gas physical conditions in the jet: the [S ii]/1673 1/6716 ratio 
increases with electron density, n e (top left); the [Nii]/[Oi] ratio mainly depends on the ionisation fraction, x e , and increases for 
increasing x e (top right); the [Sn]/Ha ratio is decreasing for increasing excitation conditions (bottom left); the [Oi]/[Sn] ratio 
depends on all the parameters but is mainly increasing with increasing temperature, T e , and density n e (bottom right). To better 
visualize the variations of the line ratios a histogram equivalent color scale is used. 



3.2. Jet physical structure: density and excitation 

The position-velocity diagrams shown in Fig.[TJsuggest a strong 
asymmetry of the gas excitation conditions in the two jet lobes. 
Low excitation lines, such as [Sn]/16731 and [Oi]/I6300, are 
much fainter in the blue lobe than in the red one, while higher 
excitation lines, such as [Nn]/16583 and Ho-/16563, are equally 



bright in the two lobes (see also the top panel of Fig.|7j). This sug- 
gests that the gas in the blue lobe is not only moving at higher 
velocity, as inferred in the previous section, but is also more ex- 
cited with respect to the gas in the red lobe. 

In order to examine the jet physical conditions we 
compute selected line ratios which are known to be good 
tracers of the gas excitation (see Fig. [5]). [S n]/16731/6716 
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Fig. 6. Position-velocity diagrams of the gas physical conditions in the jet (a histogram equivalent color scale is used). The electron 
density, n e , in units of 10 4 crrT 3 (top left), the ionisation fraction, x e (top right), the temperature, T e , in units of 10 4 K (bottom right), 
and the total density, n#, in units of 10 4 crrr 3 (bottom left) are calculated for each pixel of the spectral images. (1 pix=0.378" in 
the spatial direction and 2.15 kms -1 in the spectral direction; angular resolution: FWHM Je , c ,, w ~0.8"~2 pixels, spectral resolution: 



FWHM,„ S ,~8.8 kms 1 ~4 pixels). The errors on the estimated values of n e 
depending on the line signal-to-noise on each pixel. 



x e , T e , iih may vary between a few percent and 50% 



probes the electron density, n e , increasing up until n e 
reaches a critical density of 2.5 10 4 ctrT 3 . On the other 
hand, [N n]/l(6548+6583)/[Oi]/l(6300+6364) (hereafter 
[Nii]/[Oi]) and [Oi]/1(6300+6364)/[Sii]^(6716+6731) (here- 
after [Oi]/[Sii]) str ongly depend on the gas ion isation fraction 
and temperature (Bacciotti & Eisloffel 1999). In particu- 
lar, the [Nii]/[Oi] line ratio is increasing with x e while 
[O i]/[S n] increases with both T e and the density. Finally, the 
[Sn]^(6716+6731)/Ha^6563 line ratio (hereafter [Snj/Hor) is 
an indicator of the e xcitation level (see, e .g., the predictions 
from shock models bv Hartigan et al.l fl994). In agreement with 
iMundt et alJ (I1987L Il99lh . we found that [S n]/Ha can be as 
high as ~2Q in the red lobe and very low (~0.2-l) in the blue 
one. This confirms that the redshifted jet is generally of low 
excitation, while the gas in the blue lobe is highly excited. 



By applying the so-c alled BE technique, described in 



iBacciotti & Eisloffell (1 19991) . to the computed [Sn]i6731/6716, 
[Nn]/[Oi], and [Oi]/[Sii] ratios pixel by pixel we obtain 
position-velocity maps of the parameters characterizing the gas 
conditions, i.e. n e , x e , T e , and n# (Fig. |6}. The method relies 
on the fact that the observed forbidden lines are collisionally ex- 
cited and that, in low-excitation conditions, and provided that 
no strong sources of ionising photons are present, the ionisa- 
tion fraction of hydrogen is tightly relate d to those of nitrogen 
and ox ygen via c harge exchange pro cesses. IBacciotti & Eisloffel 
and, later, Nisini et al. (2005), verified that the lines used 



in the technique are emitted in a region of similar t emperature, 
density, and ionisation fraction. Hartigan & Morsel (120071) fur- 
ther confirmed the validity of the method by applying spectral 
diagnostic techniques to infer the physical properties of the HH 
30 jet. With this technique one can retrieve the gas physical con- 
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Fig. 7. Variations of the gas physical conditions along the jet and in the low (LVC, black diamonds) and the high (HVC, red 
diamonds) velocity components. The points are a median of the values shown in Fig. [6] over the HVC and the LVC. From top to 
bottom: (a) intensity profiles of the [Sn], [Oi], [Nn], and Ha lines integrated over their velocity profile and normalized to their 
peak intensity (b) the electron density, n e , in units of 10 3 cm 4 ; (c) the ionisation fraction, x e ; (d) the temperature, T e , in units of 10 4 
K; (e) the total density, n#, in units of 10 4 cm 4 . 



ditions comparing the observed line ratios with the ones com- 
puted for a grid of values of n e , x e , and T e . The errors affect- 
ing the estimated parameters are due mainly to the e rrors on the 
measured li ne fluxes, i.e. on the signal-to-noise (see Podi o et alj 
2006, 120091 for a discussion on the computation of the errors). 
Note that the diagnostic technique uses ratios between differ- 
ent species and, thus, the obtained x e and T e may depend on 
the chosen set of e lemental abundances. However, as stressed in 
iPodio etalJ d2006). the relative variations of the parameters (i.e. 
in different knots and/or velocity components) do not depend 
on th e choice of an abundance set (see Fig. 1 of Podi o et all 
2006) which, therefore, can be thought of as a model parame- 
ter. M oreover, the most recent det erminations by Asplun d et alj 
(120051) (Solar) and lEsteban et all d2004 (Orion) are in good 
agreement and produce diagnostic results that differ by not more 
than 15%, which in most cases is smaller than the uncertainty 
derived by measurements errors. Since the inferred physical pa- 
rameters are used also to determine the calcium gas-phase abun- 
dance with respect to S olar (see Sect. 1 4.21 . we adop t the Solar 
abund ances estimated bv lAsplund et al.l d2005l) . as in lPodio et al] 
(2009). We did not apply to the spectra a correction for extinc- 



tion, but we note that, despite the fact that the source is highly 
extincted by the surrounding optically thick disk (A y~6.5 mag 
or higher. Uones & Cohenll 19861: 1 Watson et al.ll2004l) , the Balmer 
de crement measured in the inner knots A and F (Har/H/?~6- 
7. iJones & Cohenl fl986) indicate that Ay is always <1 along 
the jet. (According to shock models Ha/H/3 may vary between 
~6-7 in low velocity shocks (V s ~20 kms 4 ) down to a mini- 
mum value of ~3 in very high velocity shoc ks (V^lOOkms" 1 ), 
lHartigan et"ai]|1994l) . On the other hand, iBacciotti & Eisloffell 
d 19991) show that, since the lines used in this diagnostic are very 
close in wavelength, errors introduced by not correcting for ex- 
tinction are at most 8%-10% for the ionization fraction and 15% 
for the temperature, as long as the visual extinction is lower 
than ~3. To check the obtained x e values we compare the ob- 
served [NiiM(6548+6583)/[Ni]/i(5198+5200) with those com- 
puted using the inferred parameters. We selected the pixels with 
high S/N since the [Ni] lines are much fainter than the lines 
used in the diagnostics. For S/N>20 we find that observed and 
computed ratios agree within 20%-30%, i.e. within the error. 

The PV diagrams in Fig.|6]show the variation of the electron 
and total density, n e and n//, the ionisation fraction, x e , and the 
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temperature, T e , both along the two jet lobes and as a function of 
the gas velocity. The PV diagrams indicate that the jet is denser 
and more excited close to the source and at high velocities in 
both lobes. The diagnostic confirms the strong asymmetry be- 
tween the two lobes, the blue one being more excited, with val- 
ues of the ionisation fraction going up to ~0.9 and temperatures 
up to 5 10 4 K. 

To have a clear view of the variation of the gas physical con- 
ditions along the two jet lobes and in the different velocity com- 
ponents we computed, at each position along the jet, the median 
of n e , x e , T e , and n# on a ~20 km s _I velocity interval centered 
on the peak velocity of the low and the high velocity compo- 
nents. The obtained values are shown in Fig. [71 while their aver- 
age over the first 5" where both the HVC and LVC are detected, 
is reported in Tab. [1] The contribution of the MVC, only barely 
visible in the blue lobe, is not significant and it is not reported in 
the figure and in the table. The error bars are the median of the 
errors on the estimated gas physical parameters over the selected 
~20 km s" 1 interval. Note that these are much larger in the blue 
lobe where the signal-to-noise is much lower. 

The electron density n e has a symmetrical trend on the first 
few arcseconds of the two lobes, decreasing from values of ~6-8 
10 3 crrT 3 close to the source to ~ 1 10 3 cm -3 at ~6" in the HVC, 
and from ~3 10 3 cirT 3 to -500 crrT 3 in the LVC. It is then 
even lower in knots C and D of the red lobe, where it reaches 
values of ~200 cm 3 . The excitation conditions, however, are 
dramatically different in the two lobes. The gas in the blue lobe 
is very hot (T e (HVC)~l-5 10 4 K, T<,(LVC)~0.6-1 10 4 K) and 
highly ionised (x c (HVC)~0.7-0.8 and x e (LVC)~0.3-0.5). On the 
contrary, in the red lobe the temperature is lower (T e ~3-5 10 3 K) 
as well as the degree of ionisation (x e ~0. 1-0.3). Finally, the total 
hydrogen density is higher in the red lobe, where it decreases 
from -2-2.5 10 4 cirT 3 in the first 1.5" down to ~5 10 3 , ~3 10 3 , 
and ~1 10 3 cirT 3 in knots A2, B2 and C, respectively; n# is 
much lower in the blue lobe, decreasing from ~1 10 4 cm~ 3 in 
knotFl to ~1 10 3 cm 3 in knot F2. 

The fact that electron and total density peak close to the 
source, where the velocity dispersion is higher, suggests that 
stronger shocks are occurring at the jet base which compress and 
heat the flow. 



4. Discussion 

In this section we derive from the inferred jet kinematical and 
physical properties crucial information on the dynamics of the 
outflows, that is the mass and momentum transported by the two 
jet lobes in their different velocity components. We then investi- 
gate the depletion of refractory species in gas phase, and finally 
we discuss the jet asymmetry and the kinematical/physical struc- 
ture of the detected velocity components comparing our results 
with those obtained for other jets from T Tauri and younger Class 
0/1 stars. 

4.1. Jet dynamics: mass and momentum flux rates 

The mass and momentum flux transported by the jet, M- ]St and 
Pjet, are fundamental parameters to understand the dynamical 
relationship between the jet and its environment, and to test the 
validity of the magneto hydro-dynamical models proposed to ex- 
plain the jet launch. For example, the models predict that the 
ratio between the rate of mass ejected into the jet (Mj et ) and ac- 
creted fro m the disk onto the star (M acc ) is fixed (Mj et / M acc ~ 
0.01 - O.l. lShu et alJl200(tlkonigl & Pudritzll2000h . 
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Fig. 8. From to p to bottom: ( a) jet radius, r jef , in arcseconds 
(as estimated by iMundt et al.l d!991l) for distances larger than 
~2" and ~3.5" in the blue and red lobe, respectively, and by 
Stanelfel dt et al.l (1 19971) for lower distances); (b, c) mass flux 
rate, Mj et , in units of 10~ 9 M yr _1 , and flux of linear momen- 
tum, Pj et , in units of 10~ 6 M yr _1 kms _1 , estimated in the high 
(HV, red points) and low (LV, black points) velocity components. 



As shown in previous studies (Hartiaan et al. 


1994; 


Bacciotti & Eisloffell 19991 iNisini et al.ll2005llPodio et all 


2006) 



an estimate of Mj et can be inferred from the determined to- 
tal density, and the measured gas velocity, Vj et , and jet ra- 
dius, rj et , as Mj e t= limn rinnr^ Vj e t, where //=1.49 is the mean 
atomic weight and m# the proton mass. The flux of linear mo- 
mentum transported by the jet, P^ t , is then easily obtained as 
/>j et = Mjet Vjet- As noted in lPodio et alj d2006l) this way of esti- 
mating the mass flux is independent on reddening and distance 
estimates. On the other hand, it assumes that the jet is uniformly 
filled at the density derived from the diagnostics within the esti- 
mated jet radius, rj et , thus providing basically an upper limit to 
Mjet- Such an effect is partially compensated for by the presence 
of regions at higher densities in the jet than those traced by the 
S nl lines, which are not taken in to account in the calculations 
Nisini et alJl2005llPodio et al.ll2006h . 

We estimate separately Mj et in the low and high velocity 
components in the two jet lobes. The absolute jet velocity, Vj et , 
is derived from the radial velocities estimated by Gaussian fit- 
ting of the [S n]/16731 line pr ofile (see Fig.[H) and corr ecting for 
the inclination angle of 65° (Eisloffel & IVl undtl[T998l) . The jet 
radius, rj et , has been estimated by IMundt et alj d!991l) for dis- 
tances larger than ~2" and ~3 / /5 in the b lue and red lobe, re- 
spectively, and bv IStapelfeldt et all d 19971) for d~071-5". The 
variation of r jeX with distance from the source, obtained by in- 
terpolating the estimates from these authors, is shown in the up- 
per panel of Fig [8] The errors on the inferred values of Mj et are 
obtained by propagating the errors affecting our estimates of the 
jet total density, n//, and of the radial velocity, V m d- One should 
keep in mind, however, that the main source of uncertainty is 
the jet radius. For example, IMundt et alj d!991l) estimated r }et 
by applying one-dimensional deconvolution method to seeing- 
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limited images of the jet in the [S n] lines and the uncertainty 
on the obtained values can be as high as 30%. The values of 
r jet obtained by IStapelfeldt etalj d 19971) through HST/WFPC2 
broadband imaging are more accurate because the jet width is 
resolved thanks to the high angular resolution offered by HST 
(~(X'l). In both cases, however, there is no spectral information, 
thus estimates of rj et relative to the different velocity compo- 
nents are not available. Here we assume that the values derived 
from the wo rks cited above are va lid f or all velocities, despite 
the fac t that iBacciotti et alJ d2000h and lLavaiiev -Fouauet et alJ 
(2000) showed that, in the case of the jet from DG Tau, the high 
velocity gas is more collimated than the low velocity component. 

The estimated values of Mj et in the HVC and in the LVC 
of the red and the blue lobe are plotted in Fig. [8] In the red 
lobe Mjet turns out to be -2-6 10~ 9 , -2-3 10~ 9 Moyr 1 in the 
high and low velocity components, respectively. In the blue lobe 
M jet (HVC)~0.2-2 10~ 8 Moyr 1 and M jet (LVC)~0.3-4 lO" 10 
M Q yr 1 . 

For what concerns the linear momentum, in the red lobe Pj et 
turns out to be -3-9 10~ 7 and -1-2 10~ 7 M yr _1 kms -1 in the 
high and low velocity components, respectively, while in the 
blue lobe it is -2-6 10~ 6 and -0.05-5 10~ 9 IVLjyr 1 kms" 1 in 
the HVC and LVC, respectively. The values of Mj et and Pj et for 
the blue lobe LVC are 2-3 orders of magnitude lower than in 
the HVC, because of the very low velocity (only a few kms -1 ) 
coupled to very low total densities (n#<5 10 3 cm -3 ), and are not 
reported in Fig. [8] In general, the much lower and higher veloc- 
ities of the LVC and the HVC in the blue lobe with respect to 
the red one, suggest less interaction with the ambient medium, 
as explained in more details in Sect. 14.31 

The average of Mj et and Pj et over the first 5" in the high 
and the low velocity components of both lobes are summarised 
in Tab. [1] The average values of Mj et indicate that, summing the 
HVC and LVC contributions, the mass outflow rate transported 
by the blue lobe (M jet (blue)~8.2 10~ 9 M yr') is only -1.3 
times higher than that transported by the red one (Mj et (red)~6.4 
10~ 9 M yr _1 ). Considering the large error on this quantity, it 
can be said that in practice the amount of mass transported by 
the jet in the two lobes is comparable, despite the large differ- 
ences in density and velocity. This situation, already found for 
the asymmetric bipol ar jet from the T Tauri star RW Aur by 
iMelnikov et alJ (l2009t) . indicates that the ejection power is sim- 
ilar on the two sides of the disk, as it is further discussed in 
Sect. 1431 

The rate of mass ejected into the jet, Mj et , cannot be com- 
pared with the mass accretion rate onto the star, M acc , because 
DG Tau B is obscured by the optically thick accretion disk, thus 
preventing veiling measurements at opt ical and/or near-infrared 
wavelengths (e.g., lHartigan et alj[l995l) . The inferred value of 
-6-8 10~ 9 M Q yr _1 is of the same order of magnitude of those 
estimated for more evolved T Tauri stars (Mj et (TTSs)~10~ 10 - 
5 10~ 8 Mpyr" 1 ) through the same method ad o pted in this 
paper (Co ffey et all [2008: M elnikov et all 12009: Maur ri et alJ 
2010 ) or also from the luminosity of optical forbidden lines 
(lHartigan etail 1 1995b . while it turns out to be lower than the 
typical values inferred for HH jets emitted by y oung Class 0/1 
sourc e s (M iet (Class 0/i)~7 10" 8 -4 1 0" 7 M yr" 1 lHartigan et all 
119941: IBacciotti & Eisloffej 119991: iPodio et al.1 120061) . The de- 
rived low values of Mj et may indicate that the bulk of the gas 
is emitting in lower excitation lines lying at far-infrared wave- 
lengths, such as the [O i] 63//m line. This hypothesis seems 
to be confirmed by prelimi nary results from the analysis of 
Herschel/PACS observations (Podi o et al.ll2010l) . 
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Fig. 9. Position-velocity diagrams of Ca gas-phase abundance 
estimated as the ratio between observed and predicted (by 
assuming Ca solar abundance) [Can]/I7291/[S n]/16731 ratios. 
There is more c alcium in gas phas e at high velocities, similarly 
to what found in lPodio et al.l (12009). and in a greement with mod - 
els predictions of dust destruction in shocks dGuillet et al.l2 009). 



The linear momentum transported by the blue lobe 
(Pjet(blue)~2.5 10~ 6 M yr _1 kms _1 ) is -3 times higher than in 
the red one (Pj et (red)~8 10~ 7 M yr _1 kms _I ), due to the multi- 
plication by the jet velocity. In order to understand the dynamical 
relationship between the jet and coaxial molecular outflows, it is 
very interesting to compare the estimated value of the momen- 
tum flux transported by the jet, with the momentum flux trans- 
porte d by the surroundi ng molecular outflow detected in CO 
lines. Mit chell et al.l (1 19941) estimated a momentum flux of -1.9 
10~ 6 M yr _1 kms _1 by assuming a distance of -160 pc and an 
inclination angle with respect to the sky of -15°. Correcting this 
value for recent more accurate estimat es of the distance (d~140 
pc) and of the inclination angle (i-25°, Eisloffel & Mundt 1998) 
one obtains / > O utflow~10~ 6 M yr _1 kms _1 . This value is compa- 
rable to the momentum estimated for the jet, indicating that the 
latter is powerful enough to drive the molecular outflow detected 
at millimeter wavelengths. How precisely the jet transfers its lin- 
ear momentum to the sur rounding medium, howe ver, remains an 
open question (see, e.g jDownes & Cabrit 20071). 

4.2. Dust reprocessing along the jet: investigating Ca, Ni, 
and Fe gas-phase abundance 

The KECK spectra of DG Tau B include also a number of emis- 
sion lines from "refractory" species, such as Fe, Ni, and Ca. 
These species are known to be highly depleted in the interstel- 
lar medium because thei r atoms may be locked onto dust grains 
(ISavage & Sembachlll996l) . Estimates of Ca, Fe, and Ni gas- 
phase abundance in Orion and £ Oph indicate that they are de- 
pleted up to fou r orders of magnitude with respect to their so- 
lar abundances (Baldwin et al ] (l991t ISavage & Sembachi ri996: 
Este ban et alj |2004). On the other hand, dust grains can be par- 
tially or completely destroyed both in the jet launching region 
and in the shocks occurring along the flow, releasing Ca, Fe, and 

Ni atoms into the gas-p hase. 

Both the models (Isella & Natta 2005) and near-infrared 
interferometric observations (Akeson et alJ 120051: lEisner et~a"D 
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120071) show that the inner part of the accretion disk around T 
Tauri and Herbig Ae/Be stars is dust-free, because dust grains are 
destroyed by the stellar radiation up to the so-called dust evapo- 
ration radius, R evp . This is located at -0.1-1.5 AU in Herbig stars 
(Eisner et al. 2007) and at -0.05-0.3 AU in T Tauri stars (Akeson 
et al. 2005). For younger Class I sources there are no direct 
estimates of R evp because of their embedded nature. However, 
we can obtain a rough estimate of the position of R evn for DG 
Tau B, assuming that the model bv llsella & Nattal d2005l) is still 
valid for such a young source. For a lumi nosity of ~0.7-0.9 L 
dJones & Cohenlfl986t IWatson et alj|2004l) and a mass accretion 
rate of at least ~10 times higher the inferred Mj et we find that 
R evp is located at about 1-6 times the position of the trunca- 
tion radius, R,, where the stellar magnetosphere truncates the 
disk. (The posi tion of R, is computed from equation (2.2a,b) of 
IShuetalJ (11994 . assuming typica l temperature a nd mass for DG 
Tau B spectral type (<K6. IWatsoiTetal] 12004). T-3000-5000 
K and M- 0.1-1 M , and stell ar magnetic field strength, B„, of 
-0.5-3 kG. lJohns-Krul l2007). Thus, if the jet originates in the 
inner region of the disk (r<R evp ) we expect Ca, Fe, and Ni to be 
present in solar abundance. On the contrary, if the jet originates 
from a region of the disk extending beyond R evp , part of the dust 
is lifted with the jet material, and these species will be observed 
in the jet to be depleted with respect to their solar abundances. 

Also the shocks occurring along the jet can affect the gas- 
phase abundances of refractory species. Theoretical models of 
C- and J-shocks showed that processes such as thermal and in- 
ertial sputtering, photoevaporation and shattering, can destroy 
from a few percent up to -90% of the dust grains depending on 
a number of parameters such as the shock velocity, the gas pre- 
shock density, the intensity of the ma gnetic field, and the size 



and structure of the dust grains (e.g ., Jones et aD 19941 ; iJonesi 



l2000llMav etd]l2000H r>raine 2003T lGuillet et alJl2009h 



Thus, investigating the gas-phase abundances of Ca, Fe, and 
Ni along jets we can test observationally the predictions of jet 
launching models and dust reprocessing in shocks. 

To this aim we analyse the spatio-kinematical distribution of 
the emission from Ca, Fe, and Ni in the jet, and we estimate the 
gas-phase abundance of calcium, [Ca/H] g£IJ (hereafter [Ca] ga! ). 
The PV diagrams in Fig. [2] show that the emission in these lines 
is confined to the high velocities both in the red and in the blue 
lobe. The velocity profiles suggest that the atoms of Ca, Fe, and 
Ni, highly depleted in the ISM, may be released in gas-phase in 
the inner disk region and/or in the shocks. In both cases higher 
abundances are expected in the high velocity component. In fact, 
on the one hand, the jet launch models predict that the gas ex- 
tracted from the inner disk region is accelerated to higher veloc- 
ities with respect to the gas extracted from the outer disk. On the 
other hand, the destruction rate in a shock is highly dependent 
on the shock velocity. Thus higher destruction rate is expected 
in the strongest shocks and at high velocities. 

We estimate the Ca gas-phase abund ance fol lowing the 
procedure il lustrated in iNisini etaD d2005l) and iPodio et all 
(2006, 120091) . that consists of comparing observed and pre- 
dicted ratios between refractory and non-refractory species, 
such as sulphur, which are expected to be present in solar 
abundance in the interstellar medium. Following these studies 
we compare observed [Can]/17291/[S ii]/16731 ratios (hereafter 
[Can]/[S n]), with those computed through the estimated 
parameters (n c , x f , and T g ) and assum ing solar abundance 
for both Ca and S dAsplundetalJ 12005). The abundance of 
calcium in gas-phase in the jet with respect to its solar value is 
given by the ratio between observed and predicted [Ca n]/[S n] 
ratios: [Ca] gflJ /[Ca] JO ; flr =([Ca n]/[S n] W([Ca n]/[S a]) predicted- 



Following Hartigan et alJ d2004l) we assume that the [Ca n] lines 
are collisionally excited and compute the level populations . 
As in previous works dNisini et al.l l2005t IPodio et al.l 2006, 
2009) we also assume that there is no calcium in the form of 
Ca°, since its ionisation potential is very low, -6.1 eV. (Note 
that emission in the Na i D lines is barely detected only at the 
peak velocity of knots Al, A3 at -6 <x). Also we assume that 
at every position along the jet, the amount of calcium ionised 
twice with respect to singly ionised calcium, i.e. Ca ++ /Ca + , is 
equal to the hydro gen ionisation fraction, x e . As explained in 
Podi o et al.l ((2009) the latter assumption is justified by the fact 
that the ionisation potential of Ca + is similar to that of hydrogen 
(~11.9eV and ~13.6eV, respectively), as are the recombination 
and collisional ionisation coefficients for temperatures lower 
than 3 10 4 K. Thus, Ca + and H° atoms can be thought of as 
undergoing the same processes when moving along the jet. 

We have no estimates of the visual extinction, thus we can- 
not correct the Ca + and S + line fluxes for reddening. This implies 
that the observed ratios (assuming Ay =0) are actually upper lim- 
its to the real values (which decrease for increasing Ay). Thus, 
the calcium depletion inferred from our analysis is a lower limit. 

Even so, the results shown in Fig. [9] show that Calcium turns 
out to be strongly depleted, its gas-phase abundance varying be- 
tween aTewrp_erc£ntand_l 5% of its solar value ([Ca/H] !O / ar -2.04 
10~ 6 . lAsplund et al]|2005l) . Interestingly, the Ca abundance is at 
a maximum close to the source and at high velocities, where the 
broad velocity profiles and high values of n f and x f suggest the 
presence of strong shocks. We also derived an estimate of the 
Ca gas-phase abundance in the blue lobe, despite the very low 
signal-to-noise of calcium emission here. It turns out that the 
observed [Can]/[S n] ratio is -1-3 times the predicted one, im- 
plying that Ca is completely released into the gas-phase in the 
blue lobe, where the very high excitation conditions (x c , T e ) in- 
dicate very high shock velocities. Simi lar results have also been 
found in the HH 111 protostellar jet dPodio et alj|2009l) and con- 
firm the dependence of elemental depletion on sho ck velocity as 
predi cted by shock dust destruction models (e.g., iGuillet et al.l 
2009). 



4.3. Investigating the origin of velocity components and 
asymmetries in stellar jets 

The presence of a low and a high velocity compo n ents in jets 
from young stars was first noted by Haman^ dl994l) . lHirth et al.l 
d 19971) analysed [Sn], [Oi], and [Nn] forbidden emission lines 
in 38 TTSs and report 58% of them showing a LVC and a 
HVC, with the LVC being always less extended than the HVC. 
Similarly, two velocity components have been detected at the 
base of HH jets emitted by Class 0/1 sources and, also in these 
cases, the LVC is less extended than the HVC. However, while 
in the fl ow from T Tauri stars the LVC is confined within -20- 
200 AU dHirthetalJI 19971) . the jets emitted by younger Class 0/1 
sources present LVC extending up to - 1000-2000 AU from the 
source dGarcia Lopez et al.l2008Ll2010l) . In agreement with these 
results the LVC from DG Tau B is observed up to -500-700 AU 

from the source; 

Moreover, iHirth et all d 19971) note that the LVC is not ob- 
served in the [N n] lines, suggesting that the LVC is less excited 
than the HVC, i.e. N is ionised only once in the low velocity gas. 
More recent deeper observations, however, have shown low ve- 
locity [Nn] emission in a few objects, such as DG Tau and CW 
Tau, allowing us to apply spectral diagnostic techniques and to 
infer the physical properties of the gas in the two velocity com- 
ponents. Therefore, we compare our results with those found for 
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DG Tau and CW Tau to investigate the properties and the origin 
of the high and the low velocity components in jets. 

The values of the physical parameters summarised in Tab. Q] 
show that in the jet from DG Tau B the HVC is denser and 
more excited than the LVC in both the blue and the red lobe. 
Higher electron densities and ion isation fractions in the HVC 
have been found a lso in DG Tau (Lavallev-Fouquetet al. 2000; 
IMaurriet al.ll2010h . and in CW Tau (ICoffev et alJl2008h . On the 
other ha nd, the analysis of the line ratios in a large sample of 
TTSs byl amann ( 1994) suggest higher temperatures but lower 
electron densities in the HVC. Combining the results from DG 
Tau B with these previous estimates, we can conclude that there 
is clear evidence that the HVC is more excited than the LVC. To 
understand if the HVC is also denser is more difficult since in 
many cases an accurate estimate of the total density is limited by 
the fact that the [S n] lines used to infer n e are saturated in one 
or both velocity co mponents, hence obtaining only lower lim- 
its to t he real density dLavallev-Fou quet et a l.l2000t ICoffev et"al] 
2008; IMaurriet alJl2010h . For iet from Class 0/1 sources, the 
electron density, estimated through near-infrared [Fe n] lines, is 
higher at high velocities along HH 1 and at low velocities in HH 
34 and HH 46-47. Since no estimate of the ionisation fraction are 
available, however, it is not possible to understand which is the 
denser component. Interestingly, higher densities and excitation 
conditions at high velocities have been noted at the position of 
the shock surfac es detected along the HH 1 1 1 jet 20-60" away 
from the source dPodio et al]l2009i) . 

The fact that the LVC shows lower excitation conditions sug- 
gests that this component is originated in an extended disk-wind 
where the gas extracted from an outer region of the disk with re- 
spect to that emitting the HVC is accelerated at lower velocities 
and give rise to slower shocks (hence, the lower ionisation frac- 
tion and temperature). This hypothesis seems to be confirmed by 
high angular resolution observations of DG Tau obtained with 
adaptive optics system and/or with the Hubble Space Telescope 
which resolve the jet perpendicular to its a xis and show that the 
HV gas is more collimated than the LV g as (Bacci otti et a U2000t 
lLavallev-Fouquet et alJl2000l;lMaurri et alj|2010t) . 

Another interesting and still unexplained characteristic of 
jets from young stars is that they are often asymmetric. 
Ikirth et alJ (119941) noted that 8 out of 15 objects for which es- 
timates of the radial velocities in both lobes are available show 
a velocity asymmetry with the ratio between the velocity of the 
fast lobe and the velocity of the slow lobe, V //V s , ranging be- 
tween 1.4 and 2.6. However, a statistical study of asymmetric 
jets, aimed to investigate the origin of such an asymmetry, is 
hampered by the fact that in most of the cases one lobe is much 
more extincted or intrinsically fainter than the other one making 
it difficult to derive the gas properties in both lobes. Our deep 
KECK/HIRES observations allows for the first time to detect 
both lobes of DG Tau B in all the emission lines which are em- 
ployed in our diagnostic techniques and thus to derive n e , x e , T e , 
and n# in both lobes. 

We compa r ed o ur results against those found by 
Meln ikov et al.1 (|2009) for the well-known asymmetric jet 
from RW Aur, finding that the jets from the two stars have very 
similar kinematical and physical characteristics. First of all, in 
both cases the blue lobes are unusually much fainter than the 
red ones. We note that in both cases the density of the blue lobe 
is lower, implying intrinsically fainter emission. The blue lobe 
is also faster in both targets, probably indicating less interaction 
with the environment on that side of the system. In more detail, 
the blue lobe is ~1.8 times faster than the red one in RW Aur 
and ~2.3 in DG Tau B, it is more ionised by a factor of ~3 in 



RW Aur and ~3.3 in DG Tau B and in both jets the fast lobe 
is less dense than the slow one by a factor ~0.3. An important 
difference between the two jets is that while the width of the 
jet from RW Aur is similar in the two lobes (rj et //rj et s ~l), 
the available images of DG Tau B indicate that the fast lobe 
is less collimated than the slow one (rj et f/rj ets ~2). The most 
interesting similarity, however, is that in both DG Tau B and RW 
Aur the mass transported by the two lobes is comparable (with 
a ratio of about 1.3 in both targets, smaller than the error on 
the estimates) despite the observed asymmetry in the physical 
quantities. Finally, the flux of linear momentum transported by 
the fast lobe is obviously larger than that transported by the 
slow one in both systems (Pj et ///jet s~l-4 for RW Aur and ~3 
for DG Tau B). We note that similar Mj et in the two jet lobes 
has also been found in other stellar jets showing asymmetries in 
the p hysical conditions as, e.g., the HH 30 jet (Bacciotti et al. 
119991) and the bipolar jet from the Herbig Ae/Be star LkHa 233 
dMelnikov et"ai]l2008l) . 

The striking similarity between the properties of these asym- 
m etric jets suggests that s imilar effects are at work. As discussed 
in iMelnikov et alJ ([2009), the fact that the global mass outflow 
rate is comparable on the two sides of the DG Tau B system ap- 
pears to indicate that, as much as in the RW Aur case, the power 
transferred by the rotating disk is the same on both sides of the 
disk, supporting the validity of a magneto-centrifugal mecha- 
nism. In this scenario, however, an asymmetric ambient medium 
can lead to an asymmetric interaction of the disk with the en- 
vironment, hence different mass load, lever arms and/or radial 
extensions of the dis k region producing the jet on the two sides 
dFerreira et alJ l2006). For example, if the ambient radiation field 
is stronger on one side, this increases the level of ionisation on 
the surface of the disk on that side, leading to enhanced mass 
load on the magnetic field lines, and/or to a larger jet launching 
region. In the case of DG Tau B, larger mass-load and interaction 
with the surrounding medium in the red lobe is suggested by the 
detec tion of a redshifted m olecular outflow in the millimeter CO 
lines (Mit chell etalil994h . 

5. Summary and Conclusions 

In this paper we investigate the kinematical and physical struc- 
ture of the HH 159 jet emitted by the young Class I source 
DG Tau B, by means of high spectral resolution observations 
acqui red with KEC K/HIRES. With respect to previous studies 
(e.g., Eisloffel & Mundt 1998) we couple the analysis of the jet 
kinematics with a detailed study of its physical and dynamical 
properties in each of its lobes and velocity components. Our 
method involves the analysis of sel ected line ratios through the 
so-called BE diagnostic technique dBacciotti & Eisloffel 1999) 
to infer the gas physical conditions (n e , x e , T e , and n//). Our main 
conclusions are as follows: 

- the velocity variations along the jet indicates that the 
ejection direction is changing due to precession and/or 
interaction with the surrounding medium. The electron 
and total density, as well as the ionisation fraction and the 
temperature, are decreasing along both lobes, due to gas 
dilution, i.e. to the jet propagation in a conical geometry, 
and/or to stronger shocks at the jet base, where the ejected 
material interacts with the dense parental cloud. 

- the lines show a complex velocity and excitation structure, 
with multiple velocity components. The high velocity 
component, HVC, is extending more than the low velocity 
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component, LVC, which is fading at 500-700 AU from the 
source, similarly to wha t previously observed for other jets 
from Class 0/1 sources dGarcia Lopez et aL 2008) a nd, o n 
smaller scales, for jets from CTTSs (e.g., iHirfh et aTf l997). 
The HVC appears to be denser and more excited suggesting 
that the LV gas is extracted from an outer region of the disk 
with respect to the HV gas and is thus accelerated at lower 
velocities. 

- the jet show a strong asymmetry: the blue lobe is faster and 
more excited than the red lobe but less dense and collimated, 
suggesting that the interaction of the ejected material with 
the ambient medium is stronger on the red-lobe side. This 
asymmetry is similar to the one observed in o ther jets, and 
in particular in RW Aur dMelnikov et al.ll2009l) . 

- despite the observed asymmetries, the mass loss rate is 
similar in the two lobes (Mj et ~6-8 10~ 9 M yr -1 ), as in other 
asymmetric jets. This means that the the power transferred 
by the rotating disk is the same on both sides of the system, 
as expected if the jet is launched by a magneto-centrifugal 
engine. The observed asymmetries can be explained in the 
framework of a MHD disk wind, if t he latter propagate s 
in an inhomogeneous ambient medium ( Ferreira et alj2 006). 

- the redshifted molecular outflow detected in the CO mil- 
limeter lines (Mitc hell et al.l ll994) supports the idea of an 
asymmetric ambient medium and a larger mass-load on the 
red-lobe side. The flux of linear momentum transported by 
the jet is comparable to that estimated for the molecular 
outflow, thus the latter can be jet-driven (Pj et ~0.8-2.5 10~ 6 
M yr _1 kms _1 ). 

- the depletion of Ca gas-phase abundance with respect to its 
solar abundance ([Ca] gfls /[Ca] so / flr ~0.01-0.15) indicates that 
most Ca atoms are locked onto dust grains. The minimum 
depletion is observed close to the source and in the HVC, 
where the strongest shocks are occurring, in agreement with 
the pr edictions from mod els of dust reprocessing in shocks 
(e.g., iGuillet et al.l F2009). The presence of dust in the jet 
may imply that the material is extracted from a region of the 
disk extending beyond the so-called dust evaporation radius, 
R eV p- However, more stringent constraints on the jet launch- 
ing region can can only be derived by analysing higher an- 
gular resolution observations, to infer Ca gas-phase abun- 
dance in the first tens of AU from the source where the gas 
has not been reprocessed by shocks, and for more unem- 
bedded sources, for which the location of R flp can be de- 
rived from near-inf rared interferometric observations (e.g., 
lAkeson et al]l2005l) . 
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